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CO ' We use thermodynamic perturbation theory to calculate the free energies and resulting phase 

diagrams of binary systems of spherical colloidal particles and interacting polymer coils in good 
solvent within an effective one-component representation of such mixtures, whereby the colloidal 
particles interact via a polymer-induced depletion potential. MC simulations are used to test the 
convergence of the high temperature expansion of the free energy. The phase diagrams calculated 
for several polymer to colloid size ratios differ considerably from the results of similar calculations 
for mixtures of colloids and ideal (non-interacting) polymers, and are in good overall agreement with 
the results of an explicit two-component representation of the same system, which includes more 
^vq . than two-body depletion forces. 
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I. INTRODUCTION 



The structure, rheology and phase behaviour of sterically stabilized colloidal dispersions are strongly affected by 
' the presence of non adsorbing polymer. Nearly fifty years ago Asakura and Oosawa realized that finite concentrations 
of polymer coils would induce an effective attraction between colloidal particles, of essentially entropic origin, the 
so-called depletion interaction pj. Since the initial colloid-polymer Hamiltonian involves only repulsive interactions 
O , between all pairs of particles, the polymer-induced effective attraction between colloids, which results from tracing 
out the polymer degrees of freedom, was referred to as "attraction through repulsion" by A. Vrij. For non- interacting 
t—I [ (ideal) polymers, the range of the depletion attraction is independent of polymer concentration, and close to the 
^ ■ polymer radius of gyration R g , while the depth of the attractive well, when two colloids touch, is proportional to 
polymer concentration. Consequently, one expects that for sufficiently high concentration, and for not too small size 
ratio q = R g /R c (where R c is the radius of the spherical colloids), the effective attraction may drive a depletion- induced 
phase separation into colloid- rich ("liquid") and colloid-poor ("gas") colloidal dispersions, similar to condensation in 
simple fluids. This phase transition, which is in fact a colloid-polymer demixing transition, was investigated by Gast 
et al. 2] , who first calculated the phase diagram from thermodynamic perturbation theory. Their findings were later 
confirmed by the Monte-Carlo (MC) simulations of Meijer and Frenkel Q, and the free volume theory of Lekkerkerker 
"j*-; ■ et al. . The predicted phase diagrams agree qualitatively with experimental findings for various colloid / polymer 
mixtures 0,0- 

More recently the question was raised of how interactions between polymer coils would affect the phase behaviour 
compared to that of ideal polymers 0, H, EJ ITol ITU ]. The early theoretical investigations into the problem were 
made at the two-component level, involving an explicit consideration of the polymer coils. However very recently 
the depletion-driven effective pair-potential between two colloidal particles in a bath of interacting polymers in good 
solvent, modelled as self-avoiding walk (SAW) polymers, was calculated by MC simulations [l^. A simple analytic 
form, with coefficients determined by the SAW polymer osmotic equation of state and surface tension, yields excellent 
agreement with the simulation data over a wide range of polymer to colloid size ratio q, and polymer concentrations 



x 

In this paper we use thermodynamic perturbation theory [l4j to calculate the phase diagram of mixtures of hard 
sphere colloids and interacting polymers within the effective one-component representation, whereby the colloidal 
particles interact via the above-mentioned depletion potential induced by the SAW polymers. The results of these 
calculations can be directly compared to the predictions of recent MC simulations of the two-component representation 
of the same system Q , which agree quantitatively with recent experiments |12| . Any discrepancies between the phase 
diagrams obtained within the effective one-component and two-component representations can then be traced back to 
more-than-two-body depletion interactions between colloidal particles, which are automatically included in the latter 
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representation, but are of course neglected in the pairwise additive effective one-component picture considered in the 
present paper. 

Thermodynamic perturbation theory was previously applied to mixtures of hard sphere colloids and ideal polymers 
within the effective one-component representation by Gast et al. 0, and improved by Dijkstra et al. [l5j . Similar 
calculations were recently published for mixtures of colloids and star polymers for several functionalities / (where / 
is the number of identical arms of the star polymer connected at the centre), again within an effective one-component 
representation, as well as within the two-component description 17). Star polymers are particularly instructive 
depletants, since upon varying the functionality, they change continuously from linear polymer (/ = 2) to hard-sphere 
like behaviour (/ — > oo). The phase diagrams obtained in ref. |l7j | for / = 2 are thus, in principle directly comparable 
to the results presented in this paper. Such a comparison will be made in section 4, but is only tentative, since 
the calculations in ref. 01 neglect the polymer concentration-dependence of the effective colloid-polymer and the 
polymer-polymer interactions which are not negligible 0, Il8l Il9j . 



II. ONE AND TWO-COMPONENT REPRESENTATIONS 



Consider a system of N c spherical colloidal particles of radius R C1 in a bath of linear polymers which are in equi- 
librium with a polymer reservoir of fixed chemical potential fj, p . The corresponding semi-grand canonical description 
is schematically represented in figure ^ The colloidal particles interact via the standard hard sphere potential, while 
each polymer is made up of L monomers or segments; segments from the same or different chains are not allowed 
to overlap. In good solvent, this excluded volume constraint is the only monomer-monomer interaction and for suf- 
ficiently large L, where chemical details become irrelevant, the interacting polymers may be accurately modelled by 
self avoiding walks on a three-dimensional lattice. The monomers, moreover, are not allowed to penetrate the hard 
sphere colloids. At finite colloid and polymer concentrations, such a binary mixture of hard spheres and interacting 
polymers poses a formidable problem to theoreticians and simulators alike. 

One coarse-graining strategy which has proved very successful is to trace out individual monomer degrees of freedom 
in the "polymers as soft colloids" paradigm, whereby the total interaction between two polymer coils, averaged over all 
monomer configurations, reduces to an effective (entropic) interaction between their centres of mass, which depends on 
polymer concentration |l8j . Similarly, one can trace out the monomer degrees of freedom, to derive a state-dependent 
effective interaction between the hard sphere colloids and the centres of mass of the polymer coils |0, 191. This 
coarse-graining procedure, which amounts to a reduction of the number of degrees of freedom of each polymer from 
0{L) to 3, leads to a two-component representation of "hard" and effective "soft" colloids, which has been exploited 
in recent MC simulations to determine the phase diagram of colloid / interacting polymer mixtures for several size 
ratios q 0- However, following the Asakura-Oosawa (AO) strategy for non-interacting polymers, one can carry the 
coarse-graining procedure one step further, by eliminating the polymer degrees of freedom altogether and determining 
the resulting depletion interactions between the colloidal particles. If this procedure is carried out in the semi-grand 
canonical ensemble, the total effective interaction energy between N c colloidal particles for any configuration r c is: 

Vcc ff (r Nc ) = V cc (r N <) -k B T In < exp[-0V cp ] >^ p ,v.t (f^) (1) 

where V cc is the direct colloid-colloid interaction energy, while V cp is the total colloid-polymer interaction; the brackets 
denote a grand-canonical average over polymer degrees of freedom at fixed polymer chemical potential, volume and 
temperature, and a given colloid configuration. In the model under consideration, V cc and V cp are the colloid-colloid 
and colloid-monomer excluded volume interactions, which are pair-wise additive. The depletion interactions between 
the N c colloidal particles are given by the second term on the r.h.s. of Eq. JQ| and are not, a priori, pair-wise additive. 
However, for sufficiently low colloid concentration, or small size ratio q, the pair-wise additive contribution dominates. 



A. Depletion potential for interacting polymers 

The depletion pair potential between an isolated pair of colloidal spheres has been determined as a function of 
the centre-to-centre distance, and over a range of interacting polymer concentrations covering the dilute and semi- 
dilute regimes, in ref. 13]. The resulting depletion pair potential is accurately reproduced by the following simple 
semi-empirical form, inspired by the Derjaguin approximation |13| : 

p Vs ( r) = l-K R cMp)Mp)(i- DTuzf ; o<*<d.(p), (2) 

otherwise 
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r is the centre-to-centre separation, while x = r — 2R C is the surface-to-surface distance between the colloids; r y w {p) 
is the polymer surface tension near a planar wall, a function of polymer bulk concentration p =_Pb determined in ref. 
[2fT | for SAW polymers; D s (p) is the range of the depletion potential given, according to ref. by: 



D s {p) 



2 lw {p) R e Jf(q) 

n(p) iz$(?-o) 



(3) 



where H(p) is the osmotic pressure of the interacting polymers taken from renormalization group calculations 21], 
namely: 



where a — 2.55, <ft — AirpR^/3, R g is the polymer radius of gyration at zero density (p — 0), and 
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The second fraction on the r.h.s of (0 accounts for finite curvature of the colloid surface, as estimated approximately 
from the AO model for ideal polymers |2(i|]: 
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which reduces to 2j^pK in the limit q — > 0. 



B. Depletion potential for ideal polymers 

For comparative purposes, we also consider the depletion pair potential induced by ideal polymers. Meijer and 
Frenkel Q showed that to a good approximation this is well represented by the Asakura-Oosawa form : 

/3Bw(r) = (-^{ 1 -|i + A(f)'} ;2R c <r<2 ( r, (6) 
[ otherwise 

where a = R c + R e /o is the radius of a sphere around the colloids from which the interpenetrable polymers are 
excluded, but with an effective R^o calculated from the insertion free energy of a single colloid in a bath of ideal 
polymer. The radius is given by eq. © 20]; for a hard wall it reduces to 2/^/n w 1.128, while it monotonically 
decreases for increasing size ratio q since the polymer can deform around the spherical colloid. For the size ratios 
considered here the curvature effects are small, on the order of a few % QjJ • 

The two pair potentials (J2J and © are always attractive. For interacting polymers the range decreases with 
increasing polymer concentration p, while the latter is constant for ideal polymers. Furthermore, at any given p, the 
depth of v s (r) is always less than that of i>id(r), so that the depletion attraction induced by interacting polymers is 
weaker than for ideal polymers. Representative examples of the depletion potentials iJSJ are shown in figure 



III. FREE ENERGY CALCULATIONS 



Our objective is to draw phase diagrams of colloid / polymer mixtures in the (r] c ,<j)p) plane, where ry c = 47rp c i?j?/3 is 
the colloid packing fraction, while 4> r p — <±np r v R? g /3 is the ratio of the polymer density in the reservoir (r) over the overlap 
density p* = 3/AirRg-, the latter conventionally separates the dilute (p p < p*) and the semi-dilute (p p > p*) regimes. 
To this end we need to calculate the free energies of the various phases. Within the semi-grand canonical ensemble 
the required free energy is related to the Helmholtz free energy F C (N C , N p , V, T) by a Legendre transformation: 

F(N C , V, T) = F{N C , p p , V, T) = F C (N C , N p , V, T) - p p N p (7) 

Since the colloidal particles interact via a hard sphere repulsion and an effective, depletion-driven attraction, the 
natural way forward is to calculate the free energies of the various phases from thermodynamic perturbation theory, 
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using the well-documented hard sphere fluid as a reference system [l4|. To second order in the high-temperature 
expansion: 

F = Fq + F\ + -F 2 

= F Q + (W N J + | [(W^) Q - (W^„)2] (8) 

where Fq = Fq(N c , V, T) is the free energy of the hard sphere fluid, (...) denotes an average over the reference system 
configurations, and Wn c is the perturbation potential energy: 

Wv„ (9) 

i<j 

with Tij — \ fi — fj\, the distance between the centres of colloids i and j; w = w(r; <jf p ) is the polymer concentration- 
dependent depletion potential J5J for interacting polymers (w = v a ) or JHJ for non- interacting polymers (w = Vn). 
We stress that within the semi-grand canonical description the depletion potential must be calculated for the polymer 
density in the reservoir, which is unequivocally determined by fixing the chemical potential fi p . 

Fq in cq. (JHJ) is calculated from the Carnahan and Starling equation of state for the hard sphere fluid [22j while 
for the hard sphere solid we adopt Hall's equation of state [23|. F% is easily expressed in terms of the hard sphere 
pair distribution function go (r) for which we adopted the Verlet-Weis parametrisation in the fluid 24] , and the form 
proposed by Kincaid and Weis for the FCC solid phase [2{|. The calculation of F% involves three and four-body 
contributions of the reference system. We have adopted the approximate expression due to Barker and Henderson 
[26|. which only involves the pair distribution function and the compressibility of the reference system. Gathering 
results: 

X = + ^ PC I ~ I {^)/ Pc J ^oM^M ■ (10) 

Note that in the solid phase, go(r) is the orientationally-averaged pair distribution function. 

In order to assess the accuracy of Fi and the convergence of the perturbation series l|10(l . we have carried out 
MC simulations to compute explicitly the fluctuation term in eq. (JSJ, which is approximated by the last term in eq. 
HI U| l. as well as the total excess free energy. The latter is most conveniently calculated by the standard A-integration 
procedure ^4|, whereby the depletion-induced perturbation Wn c is gradually switched on, resulting in: 



r (w Nc : 

Jo 



F(\ = 1) = F(\ = Q)+ I (W Nc ) x d\ (11) 

where F(X — 1) is the required free energy of the fully interacting colloid / polymer mixture, F(X = 0) = Fq is the 
free energy of the hard sphere reference system, and (Wn c )\ is the statistical average of the perturbation weighted by 
the Boltzmann factor appropriate for a system of particles interacting via the hard sphere repulsion and the partially 
switched on depletion potential Xw(r). The calculation of the free energy F hence involves several MC simulations 
to determine (Wn c )\ for a series of discrete values of A|23, typically X — nx 0.05 (1 < n < 20). 

The convergence of the perturbation series i|10|) is illustrated in figure [3] for the depletion potential J2J) and a size 
ratio q = 0.67. As expected, the convergence of the series is faster when the size ratio q is larger, the polymer 
concentration is lower and the colloid packing fraction higher. 



The accuracy of each term in the series, and of the truncated sum, were tested by MC simulations of periodic 
samples of N c — 108 colloidal particles. Representative results are presented in tables I and II. Table I shows that 
values of F\ from eq. IjlQI) are very close to the simulation data, both for the fluid and the solid. The Barker-Henderson 
approximation underestimates the absolute value of F2, by less than a factor of two in the fluid phase, and by a much 
larger factor in the FCC crystal phase, where it is totally inadequate. However, as shown in Table II, the sum of the 
first three terms of the perturbation series l|10fl yields a total free energy which is surprisingly close to the "exact" MC 
results from the A-integration. Similar comparisons show that the predictions of perturbation theory are very reliable 
both for larger and for somewhat smaller values of q (say for q > 0.3), but that the predictions rapidly deteriorate 
for small q, corresponding to narrow potential wells, as expected. This failure will be illustrated in the case q = 0.1 
at the end of the following section. 
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IV. PHASE DIAGRAMS 

Once the free energies of the fluid and FCC solid phases have been calculated from thermodynamic perturbation 
theory, as explained in the previous section, the phase diagrams can be calculated using the standard double-tangent 
construction. Since the initial two-component system is athermal, the depletion potentials @ and © are purely 
entropic, so that the temperature scales out in the Boltzmann factor, and the resulting phase diagrams are independent 
of temperature. The phase diagrams for mixtures of colloids and interacting polymers in the (j] c ,<f>p) plane are shown 
in figure 01 for four values of the size ratio q, in the range 0.1 < q < 1. T he p hase diagrams look superficially similar 
to earlier results obtained for mixtures of colloids and ideal polymers [3.ll5j| or star polymers 01 • ^ n particular, for 
the smaller size ratios, the fluid-fluid phase separation is metastable, and preempted by phase coexistence between a 
high density solid and a single low density fluid phase. Such a behaviour is a typical signature of "narrow" potential 
wells like those pictured in figure 0, 0| . For larger size ratios a stable fluid-fluid phase separation appears with a 
critical point and a triple point, and the resulting phase diagrams are not unlike those of simple atomic systems in 
the density-temperature plane (with l/<fi r p playing the role of T). 

The corresponding phase diagrams for mixtures of colloids and ideal polymers, calculated using the ideal depletion 
potential 10 are shown figureElfor comparison. While they look qualitatively similar to those for interacting polymers 
in figure^] there are a number of striking quantitative differences. Because the depletion attraction for ideal polymers 
(eq. H3) is stronger than that for interacting polymers (eq. (J5J) for the same polymer concentration, the fluid- fluid 
phase separation becomes stable at a larger (jf p for the interacting than for the non-interacting polymers. While the 
phase diagrams for q ~ 0.1 are fairly close, the differences grow with increasing q. For q ~ 1 the critical point in figure 
[S](ideal case) is at (r/ c = 0.18,0^ = 0.48) compared to (0.25,1.21) in the interacting case, while the triple points are at 
(0.47,0.92) and (0.43,1.42) respectively, indicating dramatic changes when going from ideal to interacting polymers. 
Also note that while the critical polymer concentration is practically independent of q (for q > 0.35) in the ideal case, 
it shifts to higher values as q increases in the interacting case. On the other hand, the critical colloid packing fraction 
decreases as q increases in the non- interacting case, while it is practically constant for interacting polymers. 

All these trends are similar to those reported recently in simulations of the two-component description of mixtures of 
colloids and SAW polymers Q . A detailed comparison between the present perturbation theory results for the effective 
one-component system, and the phase diagrams determined for the two-component representation is made in figure 
The agreement between the simulation data for the two-component representation and the predictions of perturbation 
theory for the effective one-component representation is seen to be reasonable, but not perfect, and to deteriorate as 
q increases. The obvious reason is that perturbation theory only includes the pair-wise additive part of the depletion 
interactions, while the two-component representation also accounts for effective many-body depletion interactions 
between colloidal particles. The fact that the phase diagrams obtained from the two-component representation are 
shifted to lower polymer concentrations relative to the prediction for the effective one-component system indicates 
that the more-than- two-body depletion interactions are overall attractive in nature |2jJ. The opposite trend was found 
in the case of ideal polymers [15j , and is consistent with the more pronounced trends in the limit of large q |l6( . 

We pointed out earlier that the convergence of thermodynamic perturbation theory is expected to deteriorate when 
the range of the attractive potential well decreases, i.e. when q decreases. To check the reliability of second order 
perturbation theory at q ~ 0.1, we have systematically computed the "exact" free energy by MC simulations, using 
the A- integration (eq. iQTJ). The phase diagrams determined with the approximate and "exact" free energies are 
compared in figure The agreement remains acceptable for the fluid-solid transition, even for q ~ 0.1, but the 
(metastable) critical point of the fluid- fluid transition is at too high a colloid packing fraction[l5|. 

Two-component simulations [jj would be very expensive for small q, because the number of polymers needed scales 
as q~ 3 . However, for sufficiently small q we don't expect many-body interactions to be important, and so our one- 
component simulation should accurately represent the true colloid / polymer system. 

A final instructive comparison is between the present results for the phase behaviour of colloid/interacting polymer 
mixtures and the results for colloid/star polymer mixtures of functionality / = 2|l7|. which reduce in fact to interacting 
linear polymers considered in the present work. The phase diagrams calculated from both depletion potentials within 
the same approximation l|10|) for the free energies are compared in figure [S] for similar size ratios o. Although, as 
explained earlier in this paper, the depletion pair potential calculated for the / = 2 star polymers [TtJ is not quite 
the same as the more accurate one we use [H , the phase-dia grams show the similar trends when compared to ideal 
polymers. 

V. CONCLUSION 

We have shown that traditional thermodynamic perturbation theory, requiring only the well-documented equations 
of state and pair distribution functions of the fluid and solid phases of the reference hard sphere system, leads to 
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reasonably accurate phase diagrams of mixtures of colloidal particles and interacting polymer coils, provided the 
appropriate concentration-dependent depletion potential between two colloidal spheres is used. As expected, the 
agreement between the predictions of the effective one-component description, and the more elaborate two-component 
description observed for low size ratios q deteriorates as q increases, due to the enhanced importance of many-body 
interaction which are neglected in the one-component picture. Nevertheless, the disagreement remains tolerable even 
at q « 1, and in view of the excellent agreement between the predictions of the two-component description and 
recent experimental data |rj|.ll2|]. we conclude that the effective one-component picture, in conjunction with standard 
thermodynamic perturbation theory, provides a reliable prediction of the phase diagrams of colloids / polymer mixtures 
in good solvent. 

A direct comparison between the phase diagrams for interacting and ideal polymers calculated at the same level 
of approximation shows considerable quantitative, and even qualitative differences between the two depletants. The 
main effect of polymer-polymer interactions is to enhance the miscibility of the colloid / polymer mixtures. Similar 
conclusions were reached by a number of different recent investigations, based on two-component approaches, inclu ding 
integral equations @, "polymers as soft colloids" 0, extensions of free- volume theory^, density functional theory |Tlj. 
and star-polymer potentials [jj} . Here we show that the differences between the two types of depletants can be 
rationalised within a one-component effective potential picture, mainly because for a given R g and density p p , the 
depletion potentials for interaction polymers are less attractive than those for interacting polymers. 

The results of the present work apply to polymers in good solvent, for which the SAW model constitutes an excellent 
representation. We plan to examine the situation where solvent quality is such that attractive forces between monomers 
can no longer be neglected |28( . Upon lowering the temperature from very high (corresponding to the SAW limit) to 
the 9 temperature, we should be able to investigate the gradual change in the phase diagrams from the fully interacting 
case to one similar to the ideal polymer limit, which have both been considered in the present paper. 
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-0.196 


0.64 


Solid 


5.578 


5.468 


5.562 


0.42 


0.22 


Fluid 


-0.316 


-0.333 


-0.339 


0.64 


Solid 


3.582 


3.368 


3.563 


0.54 


0.22 


Fluid 


-0.434 


-0.468 


-0.483 


0.64 


Solid 


1.629 


1.279 


1.605 



TABLE II: Truncated free energy densities Fq + F\ + F2, in reduced units (see table I), as ob- 
tained using equation l|10|) and MC simulation, and total free energy of the system, as obtained from 
A-integration (equation (jll|l L Whereas the results of MC simulations are accurate only in the fluid 
phase, thoses obtained using equation 110JI arc in good agreement with A-integration in both phases. 
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FIG. 1: Schematic representation of the semi-grand canonical ensemble. N c colloidal particles in 
a volume V are in osmotic equilibrium with a polymer reservoir of fixed chemical potential /i p . 
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FIG. 2: Depletion potential between two colloids as a function of the surface to surface distance, shown for 
interacting polymers (w — v s from eq. (J2J) characterised by q = 0.34,0.67 and 1.05, for <\> r v = 1.16 (a), 
and by 4> r p = 0.58,1.16 and 2.32 for q — 0.67 (b). The dashed lines indicate the ideal potential (w = 
from eq. ©) for the intermediate cases. For a given <jf and q, Vid(r) is always more attractive than v a (r). 
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FIG. 3: Ratio of the first two perturbative 
ing fraction rj c and the polymer density in 
fluid phase and the crosses (r] c > 0.5) for 



terms F2/F1, as a function of the colloid pack- 
thc reservoir <p^. The diamonds are for the 

the solid phase. The size ratio is q = 0.67. 
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FIG. 4: Phase diagrams of colloid / interacting polymer mixtures as obtained from perturbation theory for the 
effective one component system with size ratios (a) q = R g /R c = 0.0905, (b) q = 0.34, (c) q — 0.67, and (d) 
q = 1.05, in the plane of the colloid packing fraction r\ c and the polymer concentration <$> T in the reservoir. F and 
S denote the stable fluid and solid (FCC) phases. F + S and F + F denote the stable fluid-solid and (meta)stablc 
fluid-fluid coexistence regions. The solid lines denote the phase boundaries for the coexistence of stable phases, 
while the dashed lines denote the metastable fluid-fluid binodal. The critical points are indicated by asterisks. 
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FIG. 5: Phase diagrams of colloid / ideal polymer mixtures as obtained from perturbation theory for the 
effective one-component system using Vid(r) (eq. ijfjj) for the same g's as in figure 01 as functions of 
the colloid packing fraction r\ c and the polymer concentration <jf p in the reservoir. Note that the criti- 
cal points (asterisks) are always at lower (jf p than the corresponding crtitical points for interacting polymers. 
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FIG. 6: Phase diagrams of colloid / interacting polymer mixtures, as obtained from perturbation the- 
ory for the effective one-component system (lines), compared to simulations of the two-component system 
(symbols) The size ratios are (a) q = 0.34, (b) 0.67 and (c) 1.05. The critical points are in- 

dicated by crosses (two-component) and asterisks (one-component). The agreement is surprisingly good. 
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FIG. 7: Phase diagram for colloid / interacting polymer mixture, as a function of the colloid pack- 
ing fraction rj c and the polymer concentration <j) r p in the reservoir, for a small size ratio (q = 
0.0905). The solid (resp. dashed-dotted) line denotes the fluid-solid binodal obtained from per- 
turbation theory (resp. from MC simulations of the effective one-component system), while the 
dashed (resp. dotted) line denotes the metastable fluid-fluid binodal obtained by the same methods. 
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FIG. 8: Phase diagrams of colloids/polymer mixtures within three different one-component 
models: ideal polymers (dashed-dotted lines, R g — 0.67), interacting polymers (solid lines, 
R g = 0.67) and star polymers with functionality / = 2 (dashed lines, R g = 0.6). The 
critical points are indicated by asterisks, while the symbol T denotes the triple points. 



